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Abstract. In this work we present a platform for testing the device performance of 
a cavity-emitter system, using an ensemble of emitters and a tapered optical fiber. 
This method provides high-contrast spectra of the cavity modes, selective detection of 
emitters coupled to the cavity, and an estimate of the device performance in the single- 
emitter case. Using nitrogen-vacancy (NV) centers in diamond and a GaP optical 
microcavity, we are able to tune the cavity onto the NV resonance at 10 K, couple 
the cavity-coupled emission to a tapered fiber, and measure the fiber-coupled NV 
spontaneous emission decay. Theoretically we show that the fiber-coupled average 
Purcell factor is 2-3 times greater than that of free-space collection; although due to 
ensemble averaging it is still a factor of 3 less than the Purcell factor of a single, ideally 
placed center. 
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1. Introduction 

Several properties of the negatively-charged nitrogen- vacancy center (NV - ) in diamond 
make it an attractive candidate as a qubit for quantum information processing 
(QIP). First, the ground-state electron-spin decoherence time exceeds a millisecond in 
isotopically purified 12 C diamond samples [TJ. Second, an excited-state manifold [21 [3] 
provides optical access to both initialize and control the electron ground-state spin [H E] 
as well as a way to create quantum entanglement through photon interference [6]. 
Finally, hyperfme coupling of the electron spin to either N or 13 C nuclear spins [TJ [8] 
provides access to a second storage qubit needed for quantum protocols based on repeat 
until success schemes |9j [10] . 

The majority of NV~ experiments have been performed on NV centers found deep 
within the diamond substrate utilizing free-space excitation and collection optics. There 
is now significant motivation to integrate these quantum emitters into optical devices to 
build scalable measurement-based optical quantum networks. To realize these networks 
coupling NV~ centers to optical cavities will be necessary to enhance emission into the 
zero-phonon line (ZPL) [TT] and for coupling ZPL emission into waveguides for efficient 
collection. NV-cavity coupling can be obtained by coupling NV centers in diamond 
nanoparticles to non-diamond microcavities [T2| IT3| HH [15], coupling NV centers in 
single-crystal diamond to non-diamond cavities (hybrid approach) [16], [T7], and by 
fabricating microcavities out of single-crystal diamond (all-diamond approach) [18] [191 
[201 EI]- l^ue to the poor spectral properties of NV centers found in nano-diamond, 
the first approach is challenging for scalable optical networks that rely on interference 
for 'which-path' erasure. Of the single-crystal diamond approaches, the all-diamond 
approach has the advantage of a higher optical field strength at the NV site which lies 
within the cavity structure, and significant enhancement of the NV emission into the 
ZPL has already been observed in all-diamond microrings p2]. The hybrid approach, 
which relies on coupling diamond NVs to the evanescent field extending out of the non- 
diamond cavity, has the potential advantage of additional device functionality provided 
by the material properties of the waveguide layer. For example, the electro-optic 
properties of the cavity material could enable switching or cavity tuning capabilities. 

A figure of merit for the coupling between a cavity mode and dipole emitter is the 
spontaneous emission enhancement (Purcell) factor F p . For an emitter at position r in 
a monolithic cavity and on resonance with the cavity mode, 
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in which Aq is the vacuum wavelength, n c is the cavity refractive index, Q is the cavity 
mode quality factor, V is the cavity mode volume, |i£(f)| is the electric field strength 
at the dipole, and |-E(r max )| is the maximum electric field strength. There are several 
key device factors that must be optimized. The optical cavity quality factor Q must 
be sufficiently large and its mode volume V sufficiently small. The quantum emitter 
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Figure 1. (a) Schematic of tapered-fiber coupling to GaP-diamond hybrid 
microcavity. (b) Device Dl: 1.5 ^tm diameter GaP disk with thickness 250 nm. (c) 
Device D2: 1.6 /im outer diameter GaP ring with thickness 150 nm. 

must be coupled spatially to the cavity mode and spectrally on resonance. Finally, 
in addition to cavity-dipole coupling, the cavity mode must be efficiently coupled to a 
useful input/output mode (i.e. a waveguide). 

When engineering optical devices for QIP applications, device testing becomes 
prohibitively difficult when all testing must be done on single emitters. This is especially 
true in the common case where the spectral and spatial properties of the emitter are 
difficult to control. In this work we develop a test platform based on emitter ensembles 
designed to evaluate coupled emitter-cavity performance. This platform is particularly 
useful for hybrid-type devices in which the material properties of the waveguide layer 
may need to be tested before investing the resources necessary to test single-emitter 
devices. 

In this work we focus on coupling an ensemble of near-surface NV~ centers in single 
crystal diamond to GaP micro-cavities fabricated through photolithographic and dry- 
etching techniques. We show that it is possible using a tapered fiber to collect emission 
preferentially from centers that are coupled to the cavity [22| 123] . We are able to 
deterministically tune the cavity resonance several nanometers onto the NV~ ZPL and 
perform lifetime measurements using the photoluminescence (PL) collected through the 
tapered fiber. The lifetime measurements can be used to evaluate NV-device coupling. 
Although we are experimentally unable to observe a significant reduction in the ensemble 
lifetime for these particular devices, we demonstrate the feasibility of the testing method 
and theoretically show that there is a significant advantage to using the tapered fiber 
over free-space collection in ensemble measurements. 

2. Hybrid GaP-diamond cavity fabrication 

The tested devices consisted of a GaP microcavity on top of a diamond pedestal with 
an ensemble of NVs near the diamond surface (Fig. [IJi). GaP has a refractive index of 
^GaP = 3.3 that is higher than that of diamond (n-^ = 2.4), and it is transparent at 
637nm, the ZPL wavelength. Thus devices made from GaP on diamond can support 
high-Q modes that couple evanescently to the NV centers near the diamond surface. 
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Experiments were performed on a a commercially available type Ha (llO)-oriented 
diamond chip (Element 6) grown by chemical vapor deposition (CVD) with a nitrogen 
concentration of [N]<1 ppm as specified by the manufacturer. The NV ensemble was 
created by N + ion implantation at 10 keV (CORE Systems) at a dose of 2 x 10 13 cm" 2 
followed by a one hour high temperature anneal at 900 °C in a H 2 /Ar forming gas. 
To maximize the fraction of NV centers in a negatively-charged state, the sample was 
further annealed in an oxygen atmosphere at 465 °C for 2 hours |24| . In this sample there 
are three mechanisms of NV formation. First, NVs are incorporated during the CVD 
growth process. The density of these NVs is expected to be less than 10 14 cm -3 [25] 
and thus their contribution to the total collected signal can be neglected. A second 
way an NV is formed is when a vacancy created by the implantation process diffuses 
and is trapped by a native nitrogen during the annealing process. The maximum NV 
density due to this process is 2x 10 15 cm -3 based on electron irradiation studies aimed at 
maximizing the NV~ yield in similar CVD material [25]. These NVs are located within 
the first 200 nm of the surface [26] and have a maximum sheet density of ~ 4 x 10 10 cm -2 . 
Finally, an NV can form when a vacancy is trapped by implanted nitrogen. This third 
type of NV is expected to be 14 ±5 nm from the surface, as estimated using SRIM |27j. 
The estimated sheet density is 10 12 cm -2 assuming a 10% conversion rate of implanted 
ions to NV~ which is typical in our devices. Even for free-space collection the majority 
of the NV PL signal is due to the implanted nitrogen. By using a tapered fiber we are 
able to further select NVs which are coupled to the cavity. We note that if an electronic- 
grade CVD diamond substrate ([N]<1 ppb) were used only the third type of NV would 
be present. 

The starting material for the GaP microcavities was a commercially available (IQE) 
wafer consisting of a 250 nm GaP layer on an 800 nm Al0.sGa0.2P sacrificial layer grown 
by molecular beam epitaxy. For one set of devices the original 250 nm thickness was 
used while in a second set the GaP layer was thinned by dry etching to 150 nm. The 
microcavities were defined photolithographically using a mask consisting of 1.5-3 ixm 
diameter disks. The resist mask (SPR-955) was reflowed on a hotplate at 180 °C and 
was used as an etch mask in an Ar/BCls/CEj reactive ion etch (RIE) chemistry. By 
varying the optical exposure time it was found that rings could also be produced due to 
diffraction around the disk mask. 

As a further step to facilitate tapered-fiber testing, a 20 //m tall, 190 /im wide 
diamond ridge defined by a SiN mask was etched on the diamond surface using an 02/Ar 
RIE chemistry. The GaP microcavities were then transferred to the diamond ridge using 
the probabilistic "drop-down" method described in Ref. [TBj. This procedure provides 
an easy shortcut for testing single devices, but eventually this must be replaced with a 
wafer-bonding approach before integrated photonic networks can be constructed. After 
transfer the devices were etched in an oxygen RIE plasma to create a 600 nm diamond 
pedestal beneath the GaP disks. This pedestal improves the lateral confinement of the 
cavity modes, decreasing radiation loss into the diamond substrate [28J. 

Scanning electron microscope images of the two devices studied for this work are 
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Figure 2. (a)Tapered-fiber photoluminescence and transmission spectra for device 
Dl. Transmission dips indicating fiber-disk coupling coincide with PL peaks, (b) 
Estimation of the coupling rate for disk (solid lines) and ring (dashed lines) devices 
as a function of device diameter for the lowest order TE and TM modes near 637 nm. 
Two coupling geometries are calculated. In (c) the fiber is in contact with the top of 
the device. In (d) it is in contact with the side of the device. 

shown in Fig. ^p,c. Device Dl consisted of a 250 nm thick, 1.5 micron diameter GaP 
disk and device D2 consisted of a 150 nm thick, 280 nm wide, ring waveguide with a 
1.6 fim outer diameter. These devices were examples of the smallest devices made and 
had cavity resonances within tuning distance of the NV~ ZPL. 

3. Tapered-fiber coupling to high-index cavities 

A 550 nm±20 nm diameter tapered fiber was used to collect NV photoluminescence 
from the hybrid GaP-diamond cavities. The tapered fiber position was controlled by 
a piezo stage and contacted to the top or side of the GaP microcavity to collect NV 
photoluminescence as shown in Fig. [Ik Room-temperature transmission measurements 
were performed to characterize the fiber-microcavity coupling. Only in the smallest 
disk devices (sub-2 /im diameter) could transmission dips be observed as shown in 
Fig. [2£i. Transmission dips ranging from a few percent to a few tens of percent could be 
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attained. The utility of using a tapered fiber for device evaluation depends critically on 
the collection efficiency, especially in the case of single centers. We have thus estimated 
the energy decay rate from the cavity mode into the waveguide as a function of device 
diameter numerically using coupled mode theory. 

We first simulated the uncoupled lowest order TE and TM cavity modes for disk 
and ring devices using FDTD (MEEP [29]). The fiber mode was calculated using well- 
known analytic expressions for a step-index circular waveguide [30]. The fiber diameter 
was taken to be 550 nm and the fiber index rif — 1.45. In Ref. |31j . an expression for 
the coupling rate between a cavity mode and fiber is derived using power conservation. 
Using this formalism, one can show that the total decay rate of energy from the cavity 
into the fiber is given by 

2 
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where E c and Ef are the cavity and fiber electric fields, respectively, and uq = 
27rc/(637nm) is the resonant frequency. The 2D integral is taken over the fiber cross 
section Af and the ID integral is along the fiber propagation axis y with E c = at 
y = ±/. We have assumed that the cavity mode is a standing wave which decays 
equally into both forward and backward propagating waveguide modes. The cavity field 
is normalized to unit energy [32], 

fi [ H*(r) ■ H c {r)dV + f e(r)E*(r) ■ E c (r)dV = 1 
Jv Jv 

and the waveguide field is normalized to unit power, 

/ (E* f (r) x H f (r) + E f (r) x H*(r)) -ydA = ±1. 
Ja s 

H c (Hf) is the cavity (fiber) magnetic field and e is the dielectric constant. 

The results of the coupled-mode calculation are shown in Fig. |2Jd for both ring and 
disk devices for a fiber in contact with the top of the device (Fig. [2J3) and for a fiber in 
contact with the side of the device (Fig. [2ji). Experimentally there is uncertainty in the 
fiber-cavity contact position, the extent of the fiber-cavity contact, the fiber diameter, 
and the the phase of the standing-waves. Even with these uncertainties the simulations 
do provide some insight into how to achieve the best fiber-device coupling. As the device 
diameter becomes smaller it appears easiest to couple to the TE modes from the side 
of the disk. The ring device in general provides better coupling due to the thinner GaP 
layer. If we assume a reasonable cavity Q of 10,000, appreciable coupling (decay rate 
due to the taper is greater than 10% of the uncoupled decay rate) only occurs for device 
diameters less than 1.5 /im. To be able to attain critical coupling in these geometries 
the device diameter would need to be sub-micron. Such small devices are desirable for 
large Purcell factors and could be fabricated using electron-beam lithography. There is 
however a fundamental limit on the minimum device diameter d m i n before Q is limited 
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Figure 3. (a) Comparison of the NV spectra with free-space collection and with 
tapered-fiber collection for device Dl. (b) Tuning spectra for device D2. Cavity-mode 
tuning with the introduction of Xenon gas can be observed. (c,d) NV spectra collected 
through the tapered fiber for D1,D2 with the cavity tuned on and off the zero-phonon 
resonance. The cavity mode of interest is labeled with an asterisk. 



by radiation loss, Qrad(^min) — 10, 000. FTDT calculations show that for GaP-diamond 
hybrid devices this diameter can be quite small: for a 150 nm GaP layer <i min = 700 nm 
and for a 250 nm layer d m i n = 650 nm. 

4. Fiber collection of the cavity-coupled zero-phonon emission at 10K 

NV photoluminescence spectra collected in free space from the center of disk Dl and 
collected through a tapered fiber coupled to the whispering-gallery modes of the cavity 
are shown in Fig. |3^. In all PL measurements the NV centers are excited with 532 nm 
excitation. In the free space spectrum, excitation and PL collection is from the top of 
the disk through a 0.55 numerical aperture objective. The NV~ ZPL line is observed 
at 637 nm with a 1.5 nm inhomogeneous linewidth. This inhomogeneity is generally 
observed in heavily implanted samples and could be a result of strain introduced during 
the implantation process. In the tapered fiber spectrum the disk is excited at its edge 
and the excitation spot is spatially separated from the tapered-fiber contact position. 
This is to reduce scatter from uncoupled NVs into the tapered fiber. The sharp cavity 
modes can clearly be observed over the background PL with a contrast greatly exceeding 
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Figure 4. (a,c) PL spectra of coupled and uncoupled NVs in device Dl, D2. Filter 
bandwidth of the collected PL for the lifetime measurement is indicated in grey. (b,d) 
Lifetime measurements of coupled and uncoupled NVs in device D1,D2. 



what can be achieved in free-space using spatial filtering [16]. 

Due to variations in device fabrication, it is unlikely that a cavity mode will initially 
be on resonance with the 637nm zero-phonon emission. Temperature tuning is not useful 
in this system because the ZPL becomes phonon broadened at temperatures greater than 
~15 K [33J. Instead we tune the cavity by injecting xenon gas into the cryostat which 
condenses on the microcavity and alters its mode volume |34J. Microdisk resonance 
tuning up to approximately 3 nm is possible for these GaP devices which have a free 
spectral range exceeding 20 nm. It is thus necessary to fabricate many devices to ensure 
a subset of devices is within the tuning range. Fig. [3Jd depicts a series of PL spectra 
from device D2 as Xe is continuously leaked into the cryostat. An increase in the PL 
intensity as the cavity is tuned onto the ZPL resonance is observed. A comparison 
between the on-resonance and off-resonance fiber spectra for both Dl and D2 are shown 
in Figs. (3}^,d showing a marked enhancement in the ZPL collection in the resonant case. 
Quality factors of 4900 and 3500 were measured for Dl and D2 respectively. 
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5. Lifetime measurements of cavity-coupled NV centers 

The NV-cavity coupling strength can be determined from changes in the spontaneous 
emission lifetime that occur when the NV ZPL is resonant with the cavity mode. The 
decay rate for a single NV center prepared in a particular excited state i may be written 

as, 

ji = 70 + F p ,i7o,zPL , (3) 

where 7 is the total uncoupled decay rate including nonradiative channels as well 
as spontaneous emission through other modes, 7o,zpl is the zero-phonon spontaneous 
emission rate in bulk diamond, and F p i is the Purcell factor. In high-quality diamond 
samples typically 70 ~ (12ns) -1 , but in the present samples we observe 70 ~ (8 ns) -1 , 
suggesting a strong non-radiative component. We have previously measured 7o,zpl ~ 
0.03 of the total NV~ photoluminescence intensity including the phonon sidebands 
which is consistent with the literature [35]. If 7o,zpl represents a fraction 5 of the total 
uncoupled decay rate 70 then the ratio of the on-resonance lifetime r c to the uncoupled 
lifetime To would be 

Te _ 7o = 1 

t 7o + ^ P ,i7zPL l+F Pji 5' 
Using the tapered-fiber collection, we measure the SE decay curve of the cavity- 
coupled NV ensemble. This measurement, combined with a theoretical consideration 
of the effects of ensemble averaging, can be used to estimate the spontaneous emission 
enhancement of a single NV center with a particular position and orientation, as is 
discussed below. Prior to performing ensemble lifetime measurements the cavity was 
tuned to the ZPL resonance. NV emission was collected through the fiber using the 
geometry depicted in Fig. |3^. PL spectra under continuous-wave excitation of NVs 
in both the disk Dl and ring D2 devices are shown in Fig. [4^l,c. For the lifetime 
measurements, a green pulsed excitation source was used, based on supercontinuum 
generation with a Ti-Sapphire laser. The excitation source had a bandwidth of 28 nm 
centered at 520nm and a repetition rate of 4.75 MHz. The ZPL emission was filtered as 
depicted by the shaded regions in Fig. |4^,c. The SE decay curves for devices Dl and D2 
were fitted to a single exponential with decay times of 7.66 ± 0.08 ns and 8.08 ± 0.18 ns 
respectively (Fig. |4|D,d). The measurement was repeated on NV centers that were not 
coupled to the cavity mode. To collect uncoupled emission the fiber was placed across 
the center of the device and the NV centers directly under the fiber were excited. As seen 
from the PL spectra in Fig. |4^l,c, in this excitation/collection geometry the ZPL line is 
weak and broad and thus shows no features of the cavity mode. Lifetime measurements 
of uncoupled NVs yielded characteristic decay times of 7.95 ± 0.08 ns and 8.38 ± 0.18 ns 
for Dl and D2 respectively. We show below that this lifetime reduction is somewhat 
smaller than what is expected theoretically. 
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6. Simulated lifetime modification 

Next, let us discuss the lifetime modification expected theoretically for a large ensemble 
of NV centers coupled to the cavity modes of device D2. Each NV center has a number 
of characteristics that affect its excited-state dynamics, including its position relative to 
the cavity mode, the orientation of the N-V axis along one of the four (111) axes, the 
principle axes of symmetry breaking for the two excited orbital states as determined by 
local strain or electric fields, and the frequency shifts of the optical transitions, again 
determined by local fields. Nonradiative decay channels may also be present, including 
decay through the singlet states (especially from the m s — ±1 spin sublevels |36j). 

To simplify the discussion, here we neglect the dynamics of the spin sublevels, and 
consider only the two orbital excited states. Such an approximation is not suitable for 
a detailed analysis of NV dynamics, but suffices to illustrate the effect of spatial and 
orientational inhomogeneity on the radiative dynamics. Under this simplification we 
can use Eq. [3] to model the decay rate for a single NV center prepared in a particular 
excited state i, and we let 7o,zpl = (0.03)(12 ns) _1 = 0.0025 ns -1 . The Purcell factor is 
given by, 

2 

fii ■ E(r 



p _ 3 n (Xo/n c ) 3 n c 
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which is similar to Eq. [T] except that we have included the factor of n c /n^ necessary 
for the hybrid geometry and also the effect of dipole alignment and spectral detuning. 
Here, is the diamond refractive index and (ii is a unit vector parallel to the relevant 
transition dipole moment. The hybrid-cavity mode volume is given by 

v _ J d 3 fe(r)\E(r)\ 2 



e(^max)|£(r ma x)| 2 ' 



in which e(r* max ) = n 2 c is the dielectric constant at the location of the field maximum [TT] . 
The NV center has two optical transitions involving the two orbital excited states with 
dipole moments orthogonal to the (111) N-V axes. The principle axes [3TJ EE] can be 
described in terms of an angle (3\ for example for an NV oriented along [1,1, l]/\/3, 
the dipole moments are /ii = u[ cos (3 + tT 2 sin (3 and fi2 = —U\ sin (3 + ui cos (3, where 
ui = [—1,-1, 2] /a/6 and u 2 = [1,-1, 0]/v2. If the angle (3 is determined by local fields 
then we expect it to be randomly distributed, and here we assume a uniform distribution. 
It is important to note that we neglect population relaxation between orbital excited 
states. This assumption is reasonable only for measurements performed at liquid-helium 
temperatures [33] . We further assume that the two orbital excited states have the same 
unmodified decay rate. In general, we must also include the effect of detuning, and here 
Aj is the detuning between the dipole transition frequency and the central frequency 
of the cavity mode, while n = u /Q is the cavity linewidth. In the measurement, we 
detect photoluminescence through a spectral filter of 0.06 nm bandwidth centered at the 
cavity frequency, and therefore we detect emission only from a subset of NV transitions 
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Figure 5. Room-temperature PL spectra collected through the tapered fiber for D2. 
Observed resonances are in reasonable agreement with a simulated modes which include 
the effect of dispersion in the GaP [39 . A similar agreement can be found for TMo 
modes for a ring with an outer diameter of 1.96 /im. Inset: Simulated mode profile 
\E(r,z,(l) = 0)| for TE^= 21 . 

that are nearly on resonance with the cavity mode. Thus we set Aj = for the present 
calculation. 

While our model assumes a simple exponential decay for each NV transition, the 
ensemble-averaged behavior is non-exponential, but can be written as, 

POO 

J(t>0)= / dCG(C)e- (70+C70 ' ZPL)t , (5) 
Jo 

where 1(f) is the measured PL decay curve, and G(Q represents a distribution of NV 
transitions with various Purcell factors, decaying at various rates, weighted by their 
relative occurrence and by their excitation and collection efficiencies. For free-space 
collection this distribution can be written as, 

G f .,(C) « P<?> 5 > Z 3 ) «, 0) S(F p (f, a, (3) - () . (6) 

Here, a is the N-V orientation vector, and p(f, a, (3) is the NV density with respect to 
position and dipole orientation. For the spatial dependence, we approximate p(z) as a 
Gaussian distribution with respect to depth below the surface, with a mean depth of 
15 nm and variance (straggle) of 5 nm. Such a distribution corresponds to a uniform 
lOkeV nitrogen ion implantation as estimated using SRIM [27], assuming that the 
resulting NV centers incorporate only implanted nitrogen, and that the implanted 
nitrogen does not diffuse during annealing. The excitation efficiency I exc (r,a, (3) 
includes a spatial dependence that is modeled as a diffraction-limited Airy disk pattern 
corresponding to 532 nm excitation through a 0.55 numerical-aperture lens (FWHM ~ 
500 nm), centered between the inner and outer edges of the ring structure. The 
excitation efficiency is also assumed to be proportional to \p-E exc \ 2 , where the excitation 
polarization E is along [1, 0, 0]. 




Probing of diamond NV ensembles coupled to GaP microcavities 



12 




Figure 6. (a) Weighted distribution of the Purcell factor F p simulated for the 
D2 microring structure with free-space collection (red) and for collection through 
a waveguide coupled only to the cavity mode of interest (blue), (b) Simulated 
photoluminescence decay curves for unmodified NV centers with 70 = (8.4 ns)^ 1 (see 
text), for NV centers in the D2 microring structure detected with free-space collection 
(red), and for waveguide collection (blue). A single-exponential fit is also shown for 
the waveguide collection case (green, dashed) with r = 7.74 ns. Simulated curves for 
devices with Q = 35000 are also included for comparison. 



For collection through the tapered fiber, Eqj6]must be modified since, ideally, only 
light emitted into the cavity mode of interest is collected. Thus we weight by the 
emission rate into the cavity mode, which is proportional to the Purcell factor: 

Cw.g.(C) oc ^p(r,a,p)I exc (r,d,p)F p (r,d,p)5(F p (r,d,p)-(). (7) 

The lowest-radial order TE cavity field was simulated for the calculation using 
FDTD (MEEP) software [2H]- A reasonably good agreement between the experimental 
and simulated resonance frequencies could be obtained for a model device with a 1.8 pm 
outer diameter as shown in Fig. [5] Using the parameters given above, setting Q = 3500 
and considering a single standing-wave cavity mode (rather than degenerate counter- 
propagating modes) we obtain the Purcell-factor histograms shown in Fig. [6^,. For this 
mode the maximum expected Purcell factor for an NV 15 nm below the surface and 
with perfect dipole alignment is -F Pjmax = 13. However, for free-space collection the 
weighted average is only (F p ) = 1.8. For collection through the tapered fiber, because 
we preferentially detect NV centers that are better coupled to the cavity mode, the 
average increases to (F p ) = 4.2. But even with the tapered fiber it is clear that the 
random spatial and angular distributions of dipole moments reduces the mean cavity 
enhancement by a large amount. 

For a more direct comparison with experiment, PL decay curves generated from 
these distributions are shown in Fig. J6Jd. Here we have assumed 70 = (8.4 ns) -1 to match 
the behavior observed in these heavily implanted samples, and 7o,zpl — 0.0025 ns -1 , 
as discussed above. Although these decay curves do not follow a single-exponential 
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time dependence, because the average decay modification is weak compared with 70, it 
would be quite difficult experimentally to detect the deviation from single-exponential 
behavior caused by variations in F p . The effective decay lifetimes obtained by fitting 
the first 40 ns of decay to a single exponential are 8.10 ns and 7.74 ns for free-space and 
tapered-fiber collection, respectively. These are approximately equal to those expected 
from the (F p ) values given above: 8.09 ns and 7.72 ns for free-space and tapered-fiber 
collection, respectively. Fig. |6Jd also includes simulations for a device with Q = 35000. 
Q's exceeding 25000 have been observed previously in much larger hybrid GaP-diamond 
devices [16] and more work is required in determining the factors limiting the Q in the 
present devices. 

The lifetime modification shown in Fig. |4j which was at the limit of detection, 
is smaller than that predicted in this simulation, and some likely causes include the 
presence of an air gap between the GaP ring [40] and the diamond surface, or non- 
ideal NV distributions. The simulations show that collection using the tapered fiber 
provides an important advantage since it selectively filters out NV centers that are 
poorly coupled to the cavity. Nevertheless, even using the tapered fiber for collection, 
ensemble averaging reduces the effective Purcell factor by a factor of approximately 3 
compared with that expected for an ideally positioned and oriented NV center. 

7. Outlook 

There is significant motivation to integrate photonics with quantum systems for 
quantum communication, quantum information interfaces, and for quantum information 
processing. In final devices, quantum-photonic "nodes" will consist of single emitters, 
however single emitters are challenging to work with when studying and developing the 
classical photonic networks. Here we presented a method combining an ensemble of 
emitters and tapered-fiber collection which can be used for fast serial testing of a large 
number of QIP devices. 

Experimentally we showed that tapered-fiber collection provides a way to obtain 
high-contrast spectra to reliably observe cavity modes. This is because the fiber enables 
the selective detection of NVs coupled to the cavity. Compared to a free-space collection, 
the measurement more directly probes the signal that would be collected by an on-chip 
waveguide, as part of an integrated photonic network. Although not exploited in this 
work, variable coupling is also possible with the tapered fiber providing a way to optimize 
the combined cooperativity parameter that includes both F p and the waveguide-cavity 
coupling efficiency. Additionally, numerical calculations indicate that for sub-micron 
GaP-diamond hybrid devices, the cavity energy decay rate into the fiber should be able 
to exceed the uncoupled decay rate allowing for efficient photon collection. 

The use of NV ensembles provides two main advantages over single centers for device 
testing. First, the collected photoluminescence can be much greater allowing for faster 
testing. Second, device characterization does not rely on nanometer-scale placement 
of single emitters. Although the devices used to test this platform did not exhibit a 
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significant Purcell enhancement, theoretically we showed that the fiber-coupled average 
Purcell factor is 2-3 times greater than that of free-space collection. Moreover, if the 
distribution of emitters is known, the ensemble device performance can also be used to 
estimate the single-emitter device performance. 

This material is based upon work supported by the Defense Advanced Research 
Projects Agency under Award No. HR001 1-09-1-0006 and The Regents of the University 
of California. 
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